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We will present our study of the flux and spectral variability of the Crab above 100 MeV on different timescales 
ranging from days to weeks. In addition to the four main intense and day-long flares detected by AGILE 
and Fermi-LAT between Sept. 2007 and Sept. 2012, we find evidence for week-long and less intense episodes 
of enhanced gamma-ray emission that we call "waves". Statistically significant "waves" show timescales of 
1-2 weeks, and can occur by themselves or in association with shorter flares. The Sept. - Oct. 2007 gamma- 
ray enhancement episode detected by AGILE shows both "wave" and flaring behavior. We extend our 
analysis to the publicly available Fermi-LAT dataset and show that several additional "wave" episodes can 
be identified. We discuss the spectral properties of the September 2007 "wave" /flare event and show that 
the physical properties of the "waves" are intermediate between steady and flaring states. Plasma instabilities 
inducing "waves" appear to involve spatial distances I ~ 10 16 cm and enhanced magnetic fields B ~ (0.5 — 1) mG. 
Day-long flares are characterized by smaller distances and larger local magnetic fields. Typically, the deduced 
total energy associated with the "wave" phenomenon (E w ~ 10 42 erg, where E w is the kinetic energy of the 
emitting particles) is comparable with that associated to the flares, and can reach a few percent of the total 
available pulsar spindown energy. Most likely, flares and waves are the product of the same class of plasma 
instabilities that we show acting on different timescales and radiation intensities. 
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1. Introduction 



The Crab Nebula (the remnant of a Supernova ex- 
plosion witnessed by Chinese astronomers in 1054) is 
powered by a very powerful pulsar (of period P = 0.33 
ms, and spindown luminosity L s d — 5 x 10 38 ergs~ 1 ) 
(see e.g., Hester [2008]). The pulsar is energizing 
the whole system through the interaction of the par- 
ticle and wave output within the surrounding Nebula 
(of average magnetic held ~ 200/xG). The resulting 
unpulsed emission from radio to gamma rays up to 
100 MeV is interpreted as synchrotron radiation from 
at least two populations of electrons/positrons ener- 
gized by the pulsar wind and by surrounding shocks or 
plasma instabilities (e.g., Atoyan & Aharonian [1996], 
Meyer et al. [2010]). The optical and X-ray brightness 
enhancements observed in the inner Nebula, known 
as "wisps", "knots", and the "anvil" aligned with 
the pulsar "jet" (Scargle [1969]; Hester [1995], [2008]; 
Weisskopf [2000]), show flux variations on timescales 
of weeks or months. On the other hand, the average 
unpulsed emission from the Crab Nebula was always 
considered essentially stable. The surprising dis- 
covery by the AGILE satellite of variable gamma-ray 
emission from the Crab Nebula in Sept. 2010 Tavani 



ct al. 2010 2011a 



Abdo et al.|20TT Buehler et al. 2010 



and the Fermi-LAT confirmation 
started a new 



era of investigation of the Crab system. As of Sept. 
2012 we know of four major gamma-ray flares from 
the Crab Nebula detected by the AGILE Gamma- 
Ray Imaging Detector (GRID) and Fermi-LAT: (1) 
the Sept-Oct. 2007 event, (2) the Feb. 2009 event, 
(S) the Sept. 2010, and (4) the "super-flare" event 
of Apr . 2011 flBuehler et al. | [2011]; |Tavani et al. 



2011b]; |Hays et al. 1 12011|; |Striani et al. | |2011a 



Striani et al. | |2011b| ; [Buehler et al. ||2012| ) . 



In this proceeding we address the issue of the 
gamma-ray variability of the Crab Nebula on different 
timescales ranging from days to weeks. We then en- 
large the parameter space sampled by previous inves- 
tigations especially for the search of statistically sig- 
nificant enhanced emission on timescales of 1-2 weeks. 
We will present a brief overview of the current knowl- 
edge on Crab's main gamma-ray flares and the results 
of a search of 7-ray enhanced emission on timescales 
of weeks in the AGILE and Fermi-LAT database and 
a brief d iscussion of our res ults. This proceeding is 
based on Striani et al. 2013 , where the data analysis 



and discussion are presented with more details. 

2. Overview of the main gamma-ray 
flares 

Four major episodes of intense gamma-ray flaring 
from the Crab Nebula have been detected by AGILE 
and Fermi-LAT 



Abdo et al. 2011 



Buehler et al. 



20121 |Striani et al. ||2011b[ |Tavam et al. ||2011a| 



Vittorini et al. ||2011| . The definition of a "flare" 
adopted in this paper is that of a single gamma-ray 
enhancement event with a risetime T\ < 1 day and 
flux F > 700 x lO^phcnr^s- 1 above 100 MeV. 
Table 1 summarizes the flaring events that we find 
by considering the AGILE and the publicly available 
Fermi-LAT database. These events show a complex 
time structure, being composed of several sub-flares 
that we classify as Fn', with n' a progressive number. 
Fig. 1 summarizes the four major flaring episodes 
with the same temporal and flux scales. The colored 
curves are indicative of the flaring behavior that in 
most cases can be represented by an exponential fit, 
characterized by rising (n) and decay (r 2 ) timescales 
as given in Table 1. We also report in Table 1 the 
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Table I Table of the flares (F > 700 x 10 -8 phcm -2 s -1 ) of the Crab Nebula found in the AGILE and Fermi data from 
Sept. 2007. 





Name 


MJD 


n (hr) 


t 2 (hr) 


Peak Flux 


B{mG) 


7* (10 9 ) 


I (10 15 cm) 


2007 


Fi 


54381.5 


22 ±11 


10 ±5 


1000 ± 150 


1.0-2.0 


2.6-4.8 


1.2-3.6 


(AGILE) 


F 2 


54382.5 


14 ±7 


6±3 


1400 ± 200 


1.1 - 2.1 


2.3 - 4.3 


0.8 - 2.2 




F 3 


54383.7 


11 ±5 


14 ±7 


900 ± 150 


1.0-2.0 


2.6-4.8 


0.8-1.7 


2009 


F 4 


54865.8 


10 ±5 


20 ± 10 


700 ± 140 


0.7-1.3 


2.6-4.8 


0.6-1.6 


(FERMI) 


Fs 


54869.2 


10 ±5 


22 ± 11 


830 ± 90 


0.8 - 1.4 


2.6 - 4.8 


0.6 - 1.6 


2010 


F e 


55457.8 


8±4 


22 ±11 


850 ± 130 


0.7-1.3 


2.5-4.7 


0.5-1.3 


(AGILE & 


F 7 


55459.8 


6±3 


6±3 


1000 ± 100 


1.4-2.6 


2.6-4.8 


0.3-0.9 


FERMI) 


F 8 


55461.9 


19 ±10 


8±4 


750 ± 110 


0.8-1.4 


2.5-4.8 


0.9-3.1 


2011 


F 9 


55665.0 


9±5 


9±5 


1480 ± 80 


1.2 - 2.2 


2.8-5.0 


0.5-1.5 


(FERMI & 


F 10 


55667.3 


10 ±5 


24 ± 12 


2200 ± 85 


1.3-2.3 


2.7-4.9 


0.6-1.6 


AGILE) 



















The timescales n and T2 are the rise and decay timescales of the flares modelled with an exponential fit respectively. 
The characteristic length of the emitting region is deduced from the relation I — c 5 n . The Lorentz factor 7* 
characterizes the adopted model of the accelerated particle distribution function dn/dy = K/a ■ £(7 — 7*), where K is 
defined in the spherical approximation, a = 1 in the spherical case, and a < 1 for cylindrical or pancake-like volumes 
reproducing the current sheet geometry. The peak photon flux above 100 MeV is measured in units of 10 -8 phcm -2 s _1 . 
These parameters are obtained, by means of a multi-parameter fit, from the following quantities (in the observer 
frame): the position of the peak photon energy, E p oc 5y* 2 B, the peak emitted power vF oc 5 4 K/al 3 B 2 y* 2 , the rise 
time ti = l/(cS), and the cooling time T2 = 8.9 x 10 3 /[(_B/Gauss) 2 7* 2 <5], assuming 5 = 1 (see text). 

data fitting physical parameters (average local mag- 
netic field B, typical particle Lorentz factor 7*, and 
characteristic size I of the emitting region of the flar- 
ing episodes). We determine these parameters and 
their uncertainties from the time constants ri and t 2 . 
Whenever applicable (2010 and 2011 events), the AG- 
ILE and Fermi-LAT data are consistent both in flux 
and spectral properties. 



3. The September-October 2007 event 
detected by AGILE 

Fig. [2] shows the 1-day lightcurve (E > lOOMeV) 
of the Crab (pulsar plus Nebula) between Sept. 24, 
2007 and Oct. 13, 2007. with a zoom of the lightcurve 
focused on the short variability episodes in the inset. 
The x 2 analysis of this data show that there are 2 
regions of enhanced emission that are more than 5 
sigma (in a 7-day timescale) above the Crab average 
flux in the AGILE data. These regions have an aver- 
age flux of ~ 450 x 10 -8 phcm -2 s -1 and a rise and 
decay time of the order of several days. These slow 
components of enhanced 7-ray emission show features 
different than those of flares that typically have rise 
and decay times of the order of 12-24 hr, and peak 
fluxes ranging from F p ^ ~ 800 x 10 -8 phcm -2 s -1 
up to Fp^o — 2500 x 10 -8 phcm -2 s -1 (as for the 
Crab super-flare of Apr. 2011). We call these slow 
component of enhanced emission waves. We indicate 
with Wi the emission from MJD ~ 54367 to MJD ~ 



54374 and with W 2 the emission from MJD ~ 54376 
to MJD ~ 54382. The 12-hr lightcurve in the inset 
of Fig. [2] shows that the 2007 peak intensity event, 
that in our previous analysis jTavani et al. 2011a 



appeared unresolved, is actually composed of three 
different flaring components, that we indicate by Pi, 
F2, and F3. Peak fluxes, rise times and decay times 
(estimated with an exponential fit) for F\, F 2l and F3 
and for W\ and W 2 are presented in Tab. [TJ The post- 
trial probability, P pos t = 1 — (1 — p) Nt , with p the 
pre-trial probability obtained from the \ 2 test, and 
N t the number of 7-day maps (trials), is P pos t > 5er 
for both W\ and W 2 . We find that the differential 
particle energy distribution function (per unit vol- 
ume) can be described by a monoenergetic function, 
dn/dy ~ K/a ■ £(7 — 7*), where 7 is the particle 
Lorentz factor, and 7* is the monochromatic value 
of the particle energy. The electron density constant 
K is defined in the spherical approximation, with 
a = 1 in the spherical case, and a < 1 for cylin- 
drical or pancake-like volumes. A power-law distribu- 
tion and/or a relativistic Maxwellian distribution were 
also shown to be consistent with the flaring data (Ta- 
vani ct al., [2011a], Striani et al. [2011b]). We adopt 
here a monochromatic distribution that deconvolved 
with the synchrotron emissivity leads to an emitted 
spectrum practically indistinguishable from the rela- 
tivistic Maxwellian shape (e.g., Striani et al. [2011b], 
Buehler et al. [2012]). This distribution is in agree- 
ment with all available gamma-ray data (for both flar- 
ing and "wave" behavior, see below), and reflects an 
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Figure 1: Gamma-ray lightcurves above 100 MeV (12-hr time bins) from the Crab (pulsar plus Nebula) detected by 
AGILE and Fermi-LAT. From top to bottom, the Sept. - Oct. 2007 event (AGILE data), the Feb. 2009 event 
(Fermi-LAT data), the Sept. 2010 event (Fermi-LAT data), and the Apr. 2011 event (Fermi-LAT data). 



important property of the flaring Crab acceleration 
process (Tavani [2013]). With our calculations we 
find for Wi a magnetic field B = (0.8 ± 0.2) mG, a 
Lorentz factor 7* = (4 ± 1) x 10 9 , and a typical emit- 
ting length range / = (0.5 — 1.5) x 10 16 cm. For F 2 , 
we find B = (1.5 ± 0.5) mG, 7* = (3 ± 1) x 10 9 , 
and I = (1.5 ± 0.7) x 10 15 cm. In our model, 
the total number of accelerated particles producing 
the gamma-ray wave/flaring behavior is in the range 
N ~ (1 - 3) • 10 38 (Afi/47r), with ACl the solid angle 
of the 7-ray emission. It is interesting to note that 
besides the differing values of the magnetic field and 
particle densities, the typical Lorentz factor and to- 
tal number of radiating particles are similar for the 
"wave" W\ and the flare F 2 . We find that this is a 
typical behavior of the transient gamma-ray emission 
that appears to be well represented by a monochro- 
matic particle distribution function with 7* ~ (3- 
-5)xl0 9 . 



4. Search for enhanced gamma-ray 
emission in the Fermi-LAT data 



Motivated by the waves found in the AGILE data, 
we searched for a similar type of enhanced gamma-ray 
emission in the publicly available Fermi-LAT data. 
With our x 2 analysis of the Fermi data we identified 
four episodes of enhanced gamma-ray emission that 
are, in a 8-days timescale, more than 5-sigma above 
the Crab average flux. We call them W3, W4, Wq and 
W7 in Fig. [3| These episodes have apparent dura- 
tions in the range from 8 to 50 days (see the complex 
marked as W 3 — W4), and an average flux in the range 
F = (350-500) x 10- 8 phcm- 2 s -1 . The event W 7 is 
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x10" 6 




Time(MJD) 



Figure 2: Lightcurve (1-day bin) of the Sept.- Oct. 2007 Crab Nebula flare detected by AGILE. In the inset the 12-hr 
bin lightcurve around the flare. This episode is characterized by a very strong variability, with waves (black line, 
marked with a W) and flares (red line, marked with an F.) 



coincident with the enhancement episode detected 1 
by Fermi-LAT (Ojha et al., 2012). Fig. [3] shows the 
detailed lightcurves of the Fermi-LAT wave episodes 
with the largest post-trial significance. Fig.[3]shows a 
remarkable episode of "wave" enhanced emission near 
MJD = 55000. The Crab was for about 50 days above 
5fi from its standard gamma-ray flux, with an average 
flux in this period F = (340 ± 6) x 10~ 8 phcm~ 2 s" 1 . 
This event appears quite complex 2 . The total episode 
that includes /*, W 3 and W4 has a time duration of 
~ 50 days, a pre-trial probability p = 2 x 10 -25 and 
a post-trial significance a pos t > 10. For each "wave" 
that we found at ha above the Crab average emis- 
sion, we estimated the rise T\ and the decay time r 2 , 
the average flux, the peak flux, the probability of ob- 
taining the given \ 2 m the nu ll hypothesis, and the 
post-trial significance. The rise and decay timescales 
are estimated with an exponential fit. The results are 
summarized in Tab. 2. 



chrotron cooling model (see also Tavani et al. [2011a], 
Vittorini et al. [2011], Striani et al. [2011b]). Tables 
1 and 2 summarize the relevant information for the 
major "flares" and "waves". We find that for flares 
lasting 1-2 days the typical length is I ~ (1 — 2) x 10 15 
cm, the density constant in the range K/a = (2 — 
8) x 10~ 9 cm~ 3 , the typical Lorentz factor 7* = (2.5- 
— 4.5) ■ 10 9 , and the local magnetic field affecting the 
cooling phase in the range B = (1 — 2) mG. As dis- 
cussed above, the total number of radiating particles, 
in case of unbeamed (5 = 1) isotropic emission, is 
N ~ (1 - 3) • 10 38 . 

For the "wave" episodes, the typical length is I > 
10 16 cm, the density constant is in the range K/a = 
(2 - 8) x nr^cm -3 , the Lorentz factor 7* = (3 - 
5) • 10 9 , and the local magnetic field in the range B = 
(0.5 — 1) mG. The total number of particles involved 
is of the same order as in the case of flares. 



5. Constraints from a synchrotron 
cooling model 

We can deduce important physical parameters of 
the enhanced 7-ray emission by adopting a syn- 



1 Due to solar panel constraints, this event was unobservable 
by AGILE. 

2 For simplicity, we use an exponential approximation for 
the "wave" emission. Admittedly, this approximation for the 
episode of Fig. |3| c entered on MJD = 55000 is not adequate. It 
is shown in Fig. 131 for illustrative purposes only. 



6. Discussion and Conclusions 

By considering AGILE and Fermi-LAT gamma-ray 
data above 100 MeV we find that the Crab produces a 
broad variety of enhanced emission. We characterize 
this enhanced emission as short timescale (1-2 day) 
"flares" and long timescale (1 week or more) "waves". 
Given the current detection level of 1-2 day enhance- 
ments (which can be extended to longer timescales 
of order of 1-2 weeks), we cannot exclude that the 
Crab is producing an even broader variety (in flux and 
timescales) of enhanced gamma-ray emission. With 
the current sensitivities of 7-ray telescopes we can ex- 
plore an important but necessarily limited range of 
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Figure 3: The most prominent wave episodes in the Fermi-LAT data with more than 5a enhancements above the 
Crab average emission. 2-day binned gamma-ray flux values (in unit of 10 -6 ph cm -2 s _1 ) above 100 MeV as a 
function of time. We notice that the event marked as /* is intermediate between flares and waves. 



Table II Table of the waves above 5a post-trial from the Crab average emission found in the AGILE and Fermi data. 



Name 


MJD 


Duration 
(days) 


(days) 


T2 

(days) 


Average Flux 
(lQ _8 phran _2 s _1 ) 


Peak Flux 
(10 _8 phcm _2 S _1 ) 


Pre-trial 
p- value 


Post-trial 
significance 


Wi 


54368-54373 


5 


2± 1 


2 ± 1 


440 ± 40 


670 ± 200 


4.5 x 10" 8 


5.0 


w 2 


54376.5-54382.5 


6 


2.5 ±1 


2± 1 


480 ± 40 


760 ± 140 


3.0 x 10~ 9 


5.5 


W 3 


54990-55008 


18 


5 ±2.5 


10 ±5 


352 ±9 


380 ± 30 


1.0 x 10" 8 


4.6 


We, 


55988-56000 


12 


5 ±2.5 


3.5 ± 1.5 


367 ± 12 


435 ± 35 


1.8 x 10~ 12 


6.2 


W 7 


56108-56114 


6 


3 ±1.5 


3 ±1.5 


431 ± 22 


450 ± 30 


1.9 x 10" 9 


5.9 



flux and spectral variations. It is interesting to note 
that what we called "flares" and "waves" (a somewhat 
arbitrary division) share the same spectral properties. 
Given the current gamma-ray sensitivities, we could 
have detected different spectral behaviors in the en- 
ergy range 50 MeV - 10 GeV. Most likely, flares and 
waves are the product of the same class of plasma in- 
stabilities that we show acting on different timcscales 
and radiation intensities. Whether or not the insta- 



bility driver of this process is truly stochastic in flux 
and timescales will be determined by a longer moni- 
toring of the Crab Nebula. The pulsar wind outflow 
and nebular interaction conditions need to be strongly 
modified by instabilities in the relativistic flow and/or 
in the radiative properties. Plasma instabilities pos- 
sibly related to magnetic field reconnection in specific 
sites in the Nebula can be envisioned. However, ev- 
idence for magnetic field reconnection events in the 
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Crab Nebula is elusive, and no optical or X-ray emis- 
sion in coincidence with the gamma-ray flaring has 
been unambiguously detected to date (e.g., Weisskopf 
et al. [2012]). 

Both the flaring and "wave" events can be at- 
tributed to a population of accelerated electrons 
consistent with a mono- chromatic or relativistic 
Maxwellian distribution of typical energy 7* ~ (2.5- 
- 5) • 10 9 . 

We also notice that the emitted total energies that 
can be deduced for the wave (5E ltW ) and flare (SE^j) 
episodes in general satisfy the relation SE 7 , W ~ SE 7 j. 
The total gamma-ray emitted energy for the wave 
episode W\ can be estimated as 5E ltW i ~ 10 41 erg. 
Therefore, the energy associated with the wave event 
Wi, taken here as an example of Crab "wave" emis- 
sion, can reach a few percent of the total spindown 
energy. 

We conclude that the Crab "wave" events are highly 
significant and quite important from the energetic 
point of view. "Waves" typically imply regions larger 
than in the case of flares, and smaller average mag- 
netic fields. Their total emitted gamma-ray energy 
can be comparable with that associated with shorter 
flares. More observations of this fascinating phe- 
nomenon are necessary to improve our knowledge of 
the flaring Crab. 
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